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> ■ Abstract 

O 

It is shown that the low energy spectrum of mesoscopic superconductors 
coupled by Josephson interaction can be probed by two-electron tunneling 
^^ ' from a normal electrode. The Andreev reflection in the NS junction of a 

r — ' normal-superconductor-superconductor double junction (NSS transistor) pro- 

^D ■ vides a unique spectroscopic tool to probe the coherent Cooper pair tunneling 

and the energy spectrum of the Josephson (SS) junction. The ground state 
■rj" \ properties are reflected in a resonant structure of the linear conductance; ex- 

cited states with an energy as low as the Josephson coupling energy lead to a 
C^ ' threshold in the nonlinear I-V characteristic. 
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The quantum states of a mesoscopic conductor (grain) connected by weak links to particle 
reservoirs are significantly influenced by charging effects. The charging energy restricts 
tunneling of electrons between the grain and the reservoirs. At a discrete, periodic set of 
values of the electrostatic potential of the grain, which can be tuned continuously by means 
of a gate, states with different grain charges may have the same electrostatic energy. In 
the case of normal metals, many-body effects lead to a non-analytical variation of the grain 
charge with gate voltage in the vicinity of these degeneracy points [00]. However, this 
behavior develops only at very low temperatures, and has not been observed experimentally 
so far. The chances for the observation of coherent charge transfer are larger if the grain 
and particle reservoirs are in the superconducting state. Here the nature of the degeneracy 
depends on the relative magnitude of the charging energy Ec and the superconducting gap 
A. If A < Ec the degeneracy occurs between states having N and A^ + 1 electrons on the 
grain. In the opposite regime, states with the numbers N and A^ + 2 may be degenerate 
and coherent Cooper pair tunneling may become observable [Q. The crossover between the 
two regimes was studied by measurements of the charge in an electron box and of the 
switching current in a Josephson double-junction p. 

The possibility of superpositions of different charge states in a superconductor gives rise 
to a number of interesting macroscopic quantum effects |§-0. Measurements of the island 
charge at low temperature and of the critical Josephson current provide information about 
the ground state energy of the system 0], but leave inaccessible the energy spectrum and 
structure of excited states |§|-§[. On the other hand, for A > Ec, the interplay between the 
Josephson coupling and charging leads to quantization of this spectrum. At the points of 
charge degeneracy the spacing between the corresponding discrete levels is of the order of 
the Josephson coupling energy Ej, which for a weak link is much smaller than the gap A. 
This hierarchy of energies should allow one to resolve the low-energy spectrum, if a suitable 
spectroscopic tool is found. 

In this paper we show that the I-V characteristic of a double tunnel junction connecting 
a superconducting lead with a normal lead through a mesoscopic superconducting grain, 
provides the needed spectroscopic tool. In the low-bias regime, the dominant mechanism of 
transport in the NS junction between the grain and the normal electrode is due to Andreev 
reflection. Under these conditions we find that the linear conductance is sensitive to the 
mixing of different charge states in the wave functions of the Josephson junction formed by 
the grain with the superconducting electrode. In addition, the nonlinear I-V characteristic 
displays a threshold voltage that allows one to identify the energy of the excited state. 

We consider a NSS transistor with a small superconducting island coupled via tunnel 
junctions to a normal left (L) electrode and a superconducting right (R) electrode. A gate 
at potential Vg is coupled capacitively to the island, and a bias voltage ±V/2 is applied 
symmetrically to the electrodes {V > 0). The charging energy of the transistor depends 
on the number of electrons A^ on the island, and also on the number of electrons N^ — N^ 
which have passed through the transistor, 

E,,iNL, N, Nj,) =Ec{^^-N^'- ^^^-^'K v . (1) 

Here C = Cl + Cr + Cg is the total capacitance of the island, i.e., the sum of the two 
junction capacitances and the gate capacitance, and Ec = e^/2C is the charging energy 



corresponding to a single electron. The effect of the gate voltage is contained in the external 
charge Qg = CgVg + {Cl — Cr)V . All the properties of the transistor are 2e-periodic in the 
external charge. We, therefore, can restrict ourselves to the range < Q^ < 2e. 

The states of the system |A'^^, N, N^ are characterized by the number of charges in the 
electrodes and the island, Nl, Nr, and A^, and furthermore, if present, by the energies of the 
excitations in the leads and the island. We start from a reference state without quasiparticle 
excitations, which we denote by |0,0,0) with charging energy Ech{0, 0,0) = Q^/2C. Other 
states participating in the two-electron transfer (Andreev reflection in the left junction and 
pair tunneling in the Josephson junction) are |0,2,— 2), | — 2,0,2), and | — 2,2,0). Their 
energies can be found from Eq. (0). The Andreev reflection process involves an intermediate 
state I — Ifc, Ip, 0), where one electron has been transferred through the NS junction onto 
the grain. The energy of this state is Ekp = Ech{—^, 1, 0) — ^^ + ep, where ^k and ep are the 
quasiparticle energies in the normal and the superconducting grain respectively. 

The Josephson coupling, which is characterized by the energy scale Ej, mixes charge 
states differing by multiples of Cooper pairs in the two superconducting electrodes. If Qg 
takes a value near e, and the bias V is small, the states |0,0,0) and |0,2,— 2) are nearly 
degenerate and hence get mixed strongly. Similarly | — 2, 0, 2) and | —2,2, 0) are mixed. 
Near the degeneracy points the eigenstates are a superposition of two states 



7/'o = a|0,0,0)+/3|0,2,-2), 
V^i = -/?|0,0,0) + a|0,2,-2) 



(2) 



with energies 
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Here the coefficients are 



a 



I-P' 



6E. 



ch 



5EI, 



El 



(4) 



and we introduced the difference in charging energy, 5Ech = -E'c/i(0, 2, — 2) — Ech{0,0,0), 
which is 5Ech = ^Eq ( ^ — 1 ) + eV . The coefficient a is close to unity if the charging 
energy of the state |0, 0, 0) lies below that of |0, 2, —2), i.e. for 5Ech > 0, and vanishes in the 
opposite limit, while f3 shows the complementary behavior. The Josephson mixing of the 
other two states leads to the following eigenstates 



^P', =a|- 2,0,2) +/3| -2,2,0), 
^i = -/5|-2,0,2)+«|-2,2,0). 



(5) 



The coefficients a and /3 are the same as for the first pair, and the corresponding energies 
are 



-^o(i) = -^o(i) - 2eK . 



(6) 



In the above consideration we neglected the effect of charging on the Josephson couphng 
constant that leads to a relatively small enhancement @,0 of Ej. 

At low bias voltages the dominant process for charge transfer across the NS junction is 



Andreev reflection. Generalizing the expression derived in Ref. |]11| for a NSN transistor, we 
can write the amplitude for this second order tunneling process between the states V'o and 
ipQ as 

p 

\Eq — Ekp Eo — EkipJ 

In this process two electrons from the states k, | and k', | of the normal electrode tunnel 
into the grain through a junction. The latter is characterized by the tunneling Hamiltonian 
with matrix elements t^p. We suppressed the spin indices, and used the relation fp ^ = 
—Vp^l between the coefficients of the Bogoliubov transformation. The energies of the virtual 
intermediate states E^p and -Efc'p where one electron has been transferred to the island enter 
the denominators. The summation in Eq. (|^) can be performed and, for C,k,^k' -^ 0, yields 
the result 

AkM^o ^ ^o) = (^Pl^Fo{tlptl>.p), ■ (8) 

This expression involves the density of states of the island u, and an average over the 
directions of the momenta p. We introduced the function 

Fo^-- ^ 

^./A2-[E,,(-l,l,0)-Eo]2 
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A-E,h{-l,l,0) + Eo 
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The rate for the Andreev reflection process is obtained by the golden rule. After sum- 
mation over the initial states k and k' one finds for low temperatures |TT| 



Here G„ is the normal state conductance of the NS junction, Rk = h/e^ ~ 25.8A;f2 is the 
quantum resistance [GuRk "^ 1), and Ngff is the effective number of parallel channels in 
the tunnel junction [^,|I^. These parameters are conveniently absorbed in the definition of 
the Andreev conductance, Ga = G'^Rx/Neff. If the applied bias V is below some threshold 
voltage VJ;h; the only transition possible at low temperatures is the Andreev reflection be- 
tween the states ipQ and ip^. The resulting current /j-es = 2er('?/'o ~^ i^'o): due to the overlap of 
the functions a and /3, shows a pronounced structure as a function of gate charge, typical for 
resonant tunneling. This is most clearly seen in the dependence of the linear conductance 
on the gate charge 
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The result is illustrated by the lower set of curves in Fig. 1. The width of the resonance is 
characterized by Ej. On this energy scale, as long as A — Ec ^ Ej, the function Fq can be 
considered constant. 

Eq. (pup was derived under the assumption that the energy A + Ech{—i, 1,0) of the 
intermediate state lies above Eq. The resonant shape changes drastically if the supercon- 
ducting gap is lowered, such that these two energies can coincide. This occurs at two values 



(12) 



If Qg lies within the window e — 6Q*g < Qg < e + 6Q*g, the ground state is no longer composed 
of even-charge states; rather the state | — 1, Ip, 0) with one electron transferred through the 
NS junction has the lowest energy PJT3|1. In the odd state the current is low; the maximal 
conductance Gres{e), see Eq. (0), is not reached. Instead, within the abovementioned 
charge window, the conductance assumes the low value G ~ Ga^GuRkY determined by 
higher order tunneling processes through both junctions ||14|| . 

Even at A > Ec — Ej/2, Andreev reflection leads solely to transitions between the states 
"00 and ipQ only, as long as the bias voltages are below a certain threshold. This threshold 
voltage is determined by the condition Vth = {Ei — Eo)/2e. Its smallest value is Vth = Ej/2e. 
If the junction capacitances are equal, Cl = Cr, we find 



eK 



th 



4 

— J 

3 



-Eci^-1 



e 



-Ec\^-l 



1 
3" 



+ '^El 



(13) 



At higher voltages, V > V^h, Andreev reflection processes can also lead to transitions between 
the other states introduced above. In addition to the transition T{ipQ -^ i/jq) we find 
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The function Fi is defined similar to Fq, but with the energy of the initial state Eq replaced 
by El. Both Fq and Fi are approximately constant on the energy scale ii^j if A — Ec ^ Ej. 
A master equation yields the probabilities Wi, Wq for the system to be in the excited and 
ground state, respectively. I.e., 



1^0 = 1-1^1. 
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(15) 



The current then is 

+ 2e[r(7Ai ^ v^i) + r(v^i ^ i,',)]w,. (16) 

The threshold dependence of the probabihty Wi on V^ — Vfh allows one to determine the 
Josephson energy Ej. This is illustrated with the help of Fig. 1. As long as the applied bias 
eV is smaller than Ej/2, which is the smallest value of eKh, the gate voltage dependence 
of the differential conductance G = dl/dV shows the resonance (^) (two lower curves in 
Fig. 1). As soon as eV > Ej/2, the probability Wi can become nonzero, and an additional 
channel for charge transfer opens up at two values of Qg (cf. Eq. (pTsD). This leads to a 
stepwise change of the conductance (two upper curves in Fig. 1). The smallest bias voltage 
at which these jumps occur enables one to determine Ej. The magnitude 6G± of the two 
jumps depends on the applied bias. If Cl = Cr we find 

6G±{V) = Ga X 
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Another way to detect the Josephson energy is a measurement of the differential conductance 
as a function of the bias voltage V^ at a fixed gate-potential, illustrated in Fig. 2. Depending 
on the value of Qg, a jump occurs at a certain threshold Vth, cf. Eq. (0); the magnitude 
(5G(Kh)isgivenbyEq.(0). 

Spectroscopy of the Josephson coupling is only possible as long as real single electron 



tunneling processes are suppressed. Similar to the case of an NSN junction |TT|, the Andreev 
reflection can get "poisoned" once such processes become possible. Poisoning occurs at a 
certain threshold bias Vpoisom at which the energy of the odd-charge state is lower than Eq. 
If the applied bias eV exceeds eV^oison by an amount as little as the level spacing of the grain, 
the probability Wo for the system to be in an odd charge state is approximately unity and 
current will drop substantially. Therefore, for the spectroscopy of the Josephson coupling, 
the condition VJ;h < Vpoison niust be met. At resonance Q^ ~ e this condition is satisfied if 
A> Ec + Ej/2. 

The resonance structure in the /(Qg)-dependence and the threshold structure in the 
/(V^)-dependence is most pronounced if the temperature is smaller than the width of these 
features, which is given by the Josephson coupling energy Ej. In addition to intrinsic ther- 
mal fluctuations there exist fluctuations due to the external electrodynamic environment. 
They persist down to the lowest temperatures, causing voltage fluctuations proportional 
to the resistance of the external circuit. Therefore, a low impedance environment creates 
the most favorable conditions for the observation of the effects discussed here. It is also 
essential that there are no quasiparticles in the grain. This requires a sufficiently large su- 
perconducting gap A > Ec + -Ej/2, a perfect BCS density of states, and not too large bias 
voltages. In exchange for these restrictions we have a well controlled theory and unambigu- 
ous predictions for experiments. We can mention that coherent tunneling of Cooper pairs 
plays a role in a number of effects. Some examples are the gate voltage dependence of the 



critical current of SSS transistors P,pi]|, and the resonant Cooper pair tunneling at finite 

6 



bias voltage P,p!5[]. The latter effect is controlled by quasiparticle- induced dissipation or 



inffuence of the environment. The spectroscopy of coherent mixing by Andreev reflection 
offers significant advantages over these examples. It is more straightforward than inferring 
the critical current from the runaway value [§, and requires only a direct measurement of 
the I-V characteristic in a convenient regime of subgap voltages. In contrast to the experi- 
ments of Haviland et al. [|15[, the suggested method does not rely on the fluctuations due to 
the hardly controllable electrodynamic environment of the junctions. 
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FIGURES 

FIG. 1. Differential conductance G/Ga of a NSS transistor as a function of the gate charge 
Qg/e^ for various values of eV/Ec = .047, .049, .051, .053 (from bottom to top, curves off-set for 
clarity). The parameters are: Ej/Ec = 0.1, Cl = Cr, and A is large compared to Eq- The two 
lower curves {eV < Ej/2) correspond to the resonance (11), the upper curves {eV > Ej/2) show 
how threshold processes affect this resonance. 

FIG. 2. Differential conductance G/Ga as a function of bias voltage eV/Ec, for various values 
of Qg/e = .99, 1.0, 1.01 (from bottom to top, curves off-set), parameters as in Fig. 1. A jump 
occurs at F = Vth-, Eq. (14). Its magnitude is given by Eq. (18): 5G^ if Qg/e > 1 — eVth/^Ec-, and 
5G- in the opposite limit. 



